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Abstract 

The antiferromagnetic three-state Potts model on the simple-cubic lattice is 
studied using the coherent-anomaly method (CAM). The CAM analysis provides the 
estimates for the critical exponents which indicate the XY universality class, namely 
a = —0.011, (3 = 0.351, 7 = 1.309 and 5 = 4.73. This observation corroborates the 
results of the recent Monte Carlo simulations, and disagrees with the proposal of a 
new universality class. 
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1 Introduction 



The three-dimensional antiferromagnetic three-state Potts model has been investigated 
intensively during the last decade. However, our understanding of its low-temperature 
properties is still far from complete. There are several controversial results suggesting 
different universality classes of its phase transition as well as different types of the low- 
temperature ordering. In the present article, we concentrate on the critical properties. 

The existence of the finite-temperature second-order phase transition from the disor- 
dered to the low-temperature phase is generally accepted, but there has been a discussion 
about its universality class. Let us review the previous results in brief. 

Banavar et al. [Q conjectured that the model under investigation belongs to the uni- 
versality class of the three-dimensional XY (0(2)) model. Ono 0] observed a Kosterlitz- 
Thouless phase with a vanishing order parameter in his Monte Carlo study, but that 
was probably a consequence of the short simulation time. While the first Monte Carlo 
estimates of the critical exponents by Wang et al. || lay somewhere between the Ising 
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and the XY universality classes, improved measurements [[| were in a good agreement 
with the later. On the other hand, using the Monte Carlo twist method, Ueno et al. || 
obtained the critical exponents which indicate a new universality class. They also found 
some evidence of a new type of the low-temperature phase. Later, Okabe estimated 
the exponents (3 and v by Monte Carlo renormalization group approach, and obtained 
results which were again in favor of the XY class. Recently, Gottlob and Hasenbusch |]7], |8j 
obtained high-precision Monte Carlo estimates of 7 and v in a very good agreement with 
the estimates for the XY model. They also measured various critical amplitudes, again 
in full agreement with the XY universality class. Estimates for critical indices are listed 
in Table 1 together with the corresponding values calculated for the 0(2)-model || [10] . 
It seems that arguments for the XY universality class predominate. On the other hand, 



Ueno et al. [11] later argued that the low-temperature phase is an incompletely ordered 



phase, and that its nature is not compatible with the XY picture. In order to help to 
clear up this controversy, we present our results based on the coherent-anomaly method. 

We use the same technique which was applied recently to the three-dimensional Ising 
model and provided accurate critical exponents. The present treatment is nothing but 
an obvious generalization of the method described in detail in the Ref. |12|]. Therefore, 



we only briefly sketch its main points in the next section. In Section 3., we present the 
results of our analysis. Concluding remarks are given in Section 4. 



2 Mean-field solutions based on a series expansion 



The coherent-anomaly method [|T^, [TJ|, [T5[] is a general approach for investigation of critical 
phenomena (see Ref. Jl6| for a recent review). It is based on the analysis of a suitable 
set of mean-field type approximations for the given system. Here, we use the variational 
series-expansion approach |T7L [18|, |H| to generate the series of approximate solutions. 



Let us consider the simple-cubic Potts model described by the Hamiltonian 



U 



J2 6 ( s ii s i) 



<i,j> 



ida s=l 



i£b s=l 



with spin variables Sj taking on three different values, say {1, 2, 3}. The first summation in 
(0) runs over all nearest-neighbor pairs of the cubic lattice, while the second and the third 
summations correspond to the interaction with external fields. We distinguish between 
the a- and 6-sublattices in order to be able to take into account the antiferromagnetic 
order appearing below the critical point. The notation hf b = /3H% ,b is used below for the 
dimensionless external fields. 

The partition function can be rewritten as follows: 

Z = E{ s j exp(-/3H) = 

J2{s abc } n*(xyz) w ( s xyz, Sx+lyz, Sx+ly+lz, $xy+lz, S xyz +1, S x+ i yz+ i, S x+ i y+ i z+ i , S xy+ i z+ i) (2) 

where the product runs only over the triples (xyz) in which all entries are either even or 
odd, and 
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w[si, s 2 , s 3 , s 4 , s 5 , s 6 , s 7 , s 8 ) = 
exp[-/3(5(si, s 2 ) + 5(s 2 , s 3 ) + 6(s 3 , s 4 ) + <5(s 4 , s x ) + 5(s 5 , s 6 ) + <J(s 6 , s 7 ) + 
8(s 7 , 8 8 ) + S(s 8 , s 5 ) + 5(si, s 5 ) + 8(s 2 , s 6 ) + <f(s 3 , s 7 ) + <5(s 4 , s 8 ))] 
x exp[ELi K( s ( s , s 1 ) + 5(s, s 3 ) + <S(s, s 6 ) + <J(s, s 8 ))] 
x exp[ELi s 2 ) + 6{s, s 4 ) + <S(s, s B ) + <*(s, sr))] • (3) 

The sum (|2|) is nothing but the partition function of a three-state vertex model defined 
on the bcc lattice, with the vertex weights w determined by @. Our strategy is to make 
use of the gauge invariance ]2D], |2T| of the vertex models to construct a set of mean-field 
solutions of the model. We parameterize the gauge transformation in the following way: 

W a {Si, S 2 , S3, S4, S5, Sq, S7, Sg) = 
S{rJ ^i ) ri-B S 2 ) r2^3,r3- B S4,r4^S6,r6- B sr ! r 7 ^8,r-8- B ss,r s w(?"l ) r 2, T 3 , r 4 , r 5 , r 6 , r 7 , r 8 ) (4) 

Wfc(si, S 2 , S 3 , S 4 , S 5 , S 6 , S 7 , Sg) = 
SfrJ -^si,r , i^4s2,r-2^3,T-3^4s4,r4-Bs 6 ,r 6 ^4s7.r 7 -Bs8^8^5,''5 U7 ( r l' r 2; r 3) r 4, ^5; r 6; r 7 ; r 8) j (5) 

where A and -B are orthogonal matrices with their first rows parameterized as 



{A 1U A 12 , A 13 } = Jl-xl- 1, 1}/V3 + Xi{0, 1, -1}/V2 + x 2 {-2, 1, 1}/V6 



{fin, 5i2, 5i 3 } = V 1 " - 1, Ij/v 7 ^ + x 3 {0, 1, ~1}/V2 + x 4 {-2, 1, l}/^6 

(6) 

Apart from the orthogonality condition, the remaining rows are arbitrary; their parame- 
terization does not effect the calculation. Alternatively, we use the parameterization 

x\ = ta cos 4>a %2 = r A sin (ft a 

x 3 = r B cos 4> B x A = r B sin 4> B (7) 

when it is appropriate. The gauge invariance consists in the fact that the partition function 
does not change when one replaces the original weights w by the transformed weights w a 
and Wb on the sublattices (of the bcc lattice) a and b, respectively. 

Following the usual procedure of the variational series-expansion method, we generate 
a formal series expansion for the transformed vertex model described by the weights w a 
and Wb without fixing the gauge parameters {x{\. There are totally 2 x 3 8 kinds of the 
vertex weights for a general three-state model with two nonequivalent sublattices. For the 
purpose of the series expansion, we classify them into 2 x 495 classes {^°' 6 }f£o induced by 
the lattice symmetry. We fix the notation such that u>Q ,b = w a ^{l, 1, 1, 1, 1, 1, 1, 1). Then, 
we calculate the formal series expansion for the free energy in powers of oof/uog and uj\/uJq. 



(See Refs. [|T2| , f22[ for technical details.) Thus, our expansion for the logarithm of the 



partition function has the form 

L 

^ 1()P (a.','. J;, i ■ ■ ' 



t l = \ logK^) + £ / n (K- 6 K' 6 }) 



n=2 
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where {f n } are homogeneous polynomials of order n, and L denotes the maximal order 
included in the expansion. Because of the large number of variables ui, this is a huge 
formula containing thousands of terms and we cannot present it here Q. The calculation 
of the series is the limiting factor in the variational series-expansion approach. In this case 
we were able to generate the series only up to the order L = 6. Nevertheless, within the 
present formulation, the approximation JF 6 includes quite large excitation encompassing 
up to 40 original spins. (Note that the graph-counting here is rather different from that 
of the usual series expansions). 

Having calculated the formal expansion for the free energy, we can return to the 
original model described by the weights w a and Wb- Thus, Tl becomes a function of the 
temperature, external fields and of the gauge parameters {x{\. Within the variational 
series-expansion method, the gauge parameters are fixed by the stationarity conditions 

_ii = (i = 1,2,3 and 4). (9) 

OXi 

Let us describe the structure of the solutions to (|^). In the high-temperature phase and 
in zero external fields, there exists only a single solution, namely Xi = 1 (i.e. ta = 
r B — 0). New solutions characterized by finite and r# appear at the critical point. 
However, our mean-field solutions JF L exhibit a nearly perfect rotational symmetry in a 
broad region around the critical point. This means that we have always = = r and 
4>a = 4>b + 7T = 0, and, moreover, JF L (r, <fi) does not depend on the angle 0. Perturbations 
of this rotational symmetry are at least of the order r 6 , and are therefore completely 
irrelevant to the CAM analysis. This property can be shown explicitly in the lowest-order 
approximation. 

The symmetric properties of our solutions reflect the restoration of the rotational 
symmetry of the model at its critical point, and are in agreement with what was observed 
in Ref. [?[]. In fact, without this property, it would be impossible to extract any reliable 
estimates for critical exponents using the CAM. 



3 CAM analysis 

The coherent-anomaly method is based on the scaling of the so-called mean-field critical 
coefficients []I3], |14], [13]. Similarly as in Ref.|18[, we expand the approximant Tl in the 
vicinity of its critical point in order to calculate the coefficients. The only difference from 
the calculation in Ref. [IHJ is that we have four gauge parameters to take into account. 
We are interested mainly in the magnetizations and the suscetibilties 

m? = — Y Q/3 = ^ (10) 
1 dhf Xl > dHf 1 J 

Naturally, these quantities are not independent. Actually, there is only one independent 
susceptibility because we have xtt — ~ Xfef/2 (k ^ I) and xti = ~Xk?- I n the same way, 
we have effectively only one critical coefficient for the magnetizations because of the rota- 
tional symmetry mentioned above. Therefore, we omit the indices % and j corresponding 
to the three Potts states as well as the sublattice indices a and j3. 



lr The series will be available upon request from fyzikomi@savba.savba.sk 
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After some straightforward calculations, we obtain the following expressions for the 
mean-field critical coefficients of the specific heat cx, the magnetization %, the suscep- 
tibility xl and the critical magnetization m c L : 

c L = (3 c L {3*(d rrP F) 2 /d rrrr F, (11) 

m L = d hx T{-§fi c L d rr pT ld m Tf 12 , (12) 

Xl = dhxTd hx T '/drrpF (13) 
and 

m C L = d^i-GPid^F/drrrrF) 1 / 3 . (14) 

Here, the symbol d x means the derivative with respect to the gauge parameter x\ or x 2 > 
dh stands for the derivative with respect to the external field, and d r is the derivative 
with respect to the radius gauge parameter. All the derivatives in (O-O) are to be 
calculated at the critical point (3 = /3£, Xi = 1 (r = ta = tb = 0). Taking different 
gauge parameters x's and external fields h results in equivalent sets of mean-field critical 
coefficients. Namely, the critical coefficients can be rescaled so that they may be equal 
to unity in the lowest order approximation L = 0. This does not affect the subsequent 
analysis, and it turns out that one is left with only one set of coefficients for each quantity. 
This is a consequence of the model symmetry. 

Now, we can estimate the true critical exponents of the model from the CAM scaling 



formulas as 12 



,f C v a/2 /\f* _fC\\ a 

cl ~ m =(|A,ir, (15) 



t* J \ t* 



, + c x (l/2-/9)/2 /u* tcWP-V 2 

~ ^ ^lA.lf-" 2 , (16) 



t* ) \ t* 

_ l ^ ( 3V-'7^, ={{AL{r , (17) 

and 

ml - i^l (^^} =(^1)--. ,. = ~,(<i-3)/3(<i-l) (is) 




where Ax, = (t* /tl) 1 ^ 2 — (tl/t*) 1 ^ 2 with t standing for the temperature or for the inverse 
temperature; t* is the exact critical value while t c L is its L th order approximation. 

In order to extract accurate estimates of critical exponents, it is necessary to know the 
critical temperature with a high accuracy. We have used the value (3* = 0.81563 obtained 
from the high-precision Monte Carlo simulation by Gottlob and Hasenbusch [0] , and fitted 
the critical exponents to the above coherent-anomaly scaling formulas. We would like to 
stress that within the CAM approach the resulting exponents always fulfill the scaling 
relations a + 2(3 + 7 = 2 and 7 = (3(5 — 1). 

We excluded the approximants for L = 0, 2 and 3 from our analysis as usual, because 
their critical behavior is the same. Actually, they represent the same approximation which 
has a Bethe-like character; it does not take into account properly even the shortest cycle 
of the lattice. This is why the critical coefficients for < L < 3 do not follow the CAM 
scaling very well for the type of the expansion used here. On the other hand, even the 
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approximation with L = 4 takes account of quite large loops, and reflects the structure 
of the lattice much better. 

Therefore, we have fitted the critical exponents to the data obtained only from the 
approximations with L = 4, 5 and 6. The values thus obtained are shown in the bottom 
of Table 2, and the corresponding CAM plots are presented in Fig. 1. Unfortunately, 
unlike in the Ising case, we are not able to estimate the error of our results, because 
we have effectively only three data points. However, the same method works extremely 
well for the Ising model, which is in fact a special case of the present system. This is 
why we believe that our estimates would be essentially unchanged also for higher-order 
approximations. Now, let us compare our results to the previously calculated exponents 
shown in Table 1. Our estimates of 7 are in a very good agreement with the the values 
of the XY universality class. The values for a and (3 are also consistent with the XY 
universality class, though the differences between various fits indicate a relatively larger 
error. 

4 Conclusion 

We have calculated the estimates of critical exponents a, /3, 7 and S of the thre-dimensional 
antiferromagnetic Potts model. Our calculation is based on the variational-expansion 
approach combined with the coherent- anomaly method. The most difficult part of the 
computation in the present investigation was the generation of the series expansion. We 
have calculated the necessary series up to the order L = 6. This relatively short series 
provides only three points for the subsequent CAM analysis, and therefore, we cannot 
estimate the error bars for our estimates. The extension of the series seems to be rather 
infeasible even for the order L = 7. Hence, the present method may not be competitive 
with the Monte Carlo simulation for the model under investigation. However, all the recent 
estimates of critical exponents calculated directly for the Potts model come from Monte 
Carlo studies. In the present situation in which similar results by different methods are 
not available, it is important to obtain at least some results independently of the Monte 
Carlo studies, in order to confirm the universality class. The agreement of our estimates 
with the XY universality class is fairly good and the discrepancy between our results and 
the values proposed by Ueno is large. We thereby conclude that our estimates corroborate 
the XY universality class, and are clearly against the new universality class as proposed 
by Ueno [|]. 

Finally, let us make a brief remark concerning the low-temperature phase. As de- 
scribed above, our mean-field solutions exhibit approximate rotational symmetry in the 
order-parameter space. From the point of view of the mean-field critical behavior, we can 
regard them as perfectly symmetric, because the deviations from the symmetry do not 
affect the calculation of the mean-field critical coefficients. However, it is the asymme- 
try which determines the direction of the sublattice-symmetry breaking. Unfortunately, 
the deviations from the perfect symmetry are so small that the differences between the 
free energies corresponding to different solutions under the stationarity conditions (Qj are 
practically negligible compared to the expected accuracy of our approximations. This 
is why we cannot determine the type of the low-temperature ordering from the present 
calculations. It is possible that this peculiar behavior of our approximations reflects the 
nature of the low-temperature phase, namely that the order parameter of the model be- 
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comes effectively continuous. This would be in line with the restoration of the rotational 
symmetry observed by Gottlob and Hasenbusch JF]. It would also elucidate the origin of 
the XY critical indices. Clearly, the interesting properties of the low-temperature phase 
deserve further investigation. 
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Table 1. 



Author 


year 


Ref. 


a 


/3 


7 


V 


Wang 


1989 


E 






1.27(5) 


0.63(4) 


Wang 


1990 


i 






1.31(3) 


0.66(3) 


Ueno 


1989 





0.15 


0.34(2) 


1.10(2) 


0.58(1) 


Okabe 


1992 


§ 




0.33(2) 




0.66(2) 


Gottlob 


1994 


E 






1.310(9) 


0.664(4) 


Gottlob 


1994 


§ 








0.663(4) 


LeGuillou 


1980 




-0.007(6) 


0.345(2) 


1.3160(25) 


0.669(2) 


LeGuillou 


1985 







0.3485(35) 


1.315(7) 


0.671(5) 



Table 1. Recent estimates for critical exponents of the antiferromagnetic three-state 
Potts model in three dimensions. We have also included the exponents calculated for 
the 3D XY model for comparison. The values in Refs. and 
estimates and the results from the e-expansion, respectively. 



10 1 are the standard RG 



Table 2. 



L 




CL 


m L 


XL 


m L 


0-3 


1.3943 


1 n 

J. . U 


1 n 

J. .U 


i n 

J. .yj 


1 n 

J. . U 


4 


1.3014 


1.09957 


1.15608 


1.21551 


1.17556 


5 


1.2998 


1.11207 


1.16766 


1.22603 


1.18680 


6 


1.2735 


1.09982 


1.23997 


1.39797 


1.29054 


data used 


a 


P 


7 


5 


5,6: 




-0.025 


0.359 


1.306 


4.63 


4,6: 




0.001 


0.344 


1.311 


4.81 


4,5,6 




-0.011 


0.351 


1.309 


4.73 



Table 2. Critical temperatures and critical mean- field coefficients as calculated for various 
approximation orders L. The critical exponents were fitted to the data from approxima- 
tions with L > 3 and the resulting estimates are shown in the bottom row of the table. 
The numbers in the first column correspond to the approximations used. 
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Figure caption 

Fig. 1. CAM scaling of the critical mean-field coefficients Ql for the specific heat, cl 
(O), magnetization, m L (+), susceptibility, xl ( n ) and for the critical magnetization, rh c L 
(x). The distance from the true critical point is measured in A L = {(3* / fli) 1 / 2 — (01/ (3*) 1 / 2 . 
Critical coefficients were rescaled so that Qo=l (see the text). 
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